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ABSTRACT 


! 

\ 

An  experimental  and  theoretical  investigation  was  conducted  to  study  the  wick  combustion  of 
Li  and  SF<s.  A  single-line  laser  induced  fluorescence  thermometry  with  Li2  as  fluorescence  species 
was  developed.  The  calibration  experiments  yielded  promising  results;  refinements,  however,  are 
needed  before  the  technique  can  be  applied  to  combustion  flame  measurements.  Wick  combustion 
of  ethanol  and  hexane  in  air  as  well  as  Li  and  SF6  was  conducted  in  a  vacuum  chamber  at 
atmospheric  and  sub-atmospheric  pressures.  The  ethanol  and  hexane  wick  diffusion  flames 
showed  flame  oscillations,  similar  to  flame  flickers  in  buoyant  jet  diffusion  flames.  The  oscillation 
frequency  was  estimated  in  the  range  5  to  8  Hz  at  the  pressures  examined.  The  flame  stand-off 
distance  was  found  to  increase  when  the  system  pressure  was  decreased.  The  wick  combustion  of 
Li  and  SF6  resulted  in  a  bright  pinkish  flame.  The  luminous  zone  appeared  quite  close  to  the  wick 
surface  at  the  condition  examined.  Near  the  "complete  consumption"  of  lithium,  the  wick  was 
burned  out  due  to  loss  of  liquid  lithium  as  a  heat  sink.  On  the  analysis,  a  conserved  scalar 
approach  was  employed  to  model  the  wick  flame.  A  single  equation  was  obtain  to  describe  the 
interface  condition  of  the  wick  combustion.  Numerical  solutions  were  obtained  for  laminar  wick 
diffusion  flames.  The  prediction  yields  similarity  profiles  for  both  ethanol-air  and  Li-SF6  wick 
diffusion  flames  although  non-similar  governing  equations  were  retained  in  the  formulation.  The 
flame  stand-off  distance  was  found  to  increase  as  the  system  pressure  was  reduced,  in  agreement 
with  experimental  observation.  The  effects  on  mass  burning  idles  due  to  fuel  composition 
variations  were  studied  employing  fuel-rich  and  product-rich  liquids  at  0.0035  MPa  in  the 
analysis.  It  was  found  that  the  product-rich  fuel  yields  a  much  lower  flame  temperature;  the  mass 
burning  rate,  however,  was  higher  than  the  pure  lithium  at  0.01  MPa.  This  is  a  result  of  the 
decrease  in  flame  stand-off  distance,  yielding  an  increase  in  the  heat  transfer  to  the  wick  surface 
when  a  product-rich  fuel  was  used.  It  is  also  noted  that  the  product  rich  liquid  has  a  lower  specific 
enthalpy  of  gasification  compared  to  fuel  rich  liquids  and  pure  lithium.  The  findings  on  fuel  mass 
burning  rates  are  interesting  -  suggesting  that  the  wick  configuration  may  achieve  a  higher  fuel 
utilization  than  the  submerged  jet.  The  future  work  for  the  next  report  period  was  also  discussed. 
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solid  angle 
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I.  INTRODUCTION 
1 . 1  General  Statements  of  the  Problem 

A  novel  liquid-metal  combustion  system  involving  the  reaction  of  lithium  (Li)  and  sulfur 
hexafluoride  (SFg)  is  studied.  The  system  is  unique  in  that  it  has  a  high  energy  density  and  a  high 
specific  energy  and  that  it  yields  condensed  phase  products  under  normal  operation  conditions. 
The  high  energy  density  and  high  specific  energy  satisfy  the  volume  and  weight  requirements  of 
propulsion  applications  and  the  condensed  phase  product  provides  a  means  of  closed  system 
operation  as  desired  by  deep-sea  propulsion  vessels;  thus  the  liquid-metal  combustion  system  is 
ideal  for  undersea  applications. 

Two  major  system  configurations  have  been  identified,  i.e.  the  submerged  jet  and  the  reactive 
heat  pipe  (or  wick-type),  e.g.  see  Hughes  et  al.  (1983).  The  submerged  jet  was  identified  as  a 
practical  energy  source  for  submersible  vehicles  (Parnell,  1987  and  1989);  however,  the  material 
compatibility  of  lithium  and  the  high  temperature  and  condensed  phase  products  of  Li  and  SF6 
combustion  makes  detailed  measurements  of  the  flame  structure  extremely  difficult  in  submerged 
jet  combustors.  The  wick  configuration,  on  the  other  hand,  provides  an  optical  access  to  study  the 
combustion  of  Li  and  SFg.  The  wick  configuration  also  has  intrinsic  merits  as  an  energy  source 
for  surface  and  under-sea  propulsion  systems. 

When  a  wick  was  employed,  liquid  Li  was  supplied  through  the  wick  action  and  heat  transfer 
was  accomplished  by  evaporation  and  condensation  of  lithium  (Hughes  et  al.,  1983  and  Faeth  et 
al.,  1978).  The  wick  system  was  shown  to  operate  at  a  lower  power  level  but  with  a  longer  time 
duration  compared  to  submerged  jet  systems.  A  submerged  jet  combustor  requires  direct  injection 
of  gaseous  SFg  into  the  liquid  bath  (Li)  at  under-expanded  conditions.  As  a  result,  the  injector  was 
in  contact  with  hot  combustion  products/lithium  and  the  injector  clogging  and  erosion  becomes  a 
problem  which  needs  to  be  addressed.  WTen  the  injector  clogged,  the  oxidant  flow  was  stopped 
and  the  combustion  was  terminated.  Restart  of  a  submerged-jet  combustor,  unfortunately,  is  very 
d  Ticult.  The  wick  combustor,  on  the  other  hand,  can  be  designed  free  from  injector  clogging. 
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This  was  accomplished  by  confining  the  high  temperature  zones  to  a  region  away  from  the  injector 
(Faeth  et  al.,  1978).  It  was  also  noted  (e.g.,  see  You  and  Faeth,  1977)  that  the  wick  system  can  be 
designed  to  have  a  fast  response  to  the  load  change,  to  be  capable  of  system  restart  and  to  have 
flexible  heat  exchanger  layout  through  the  heat  pipe  operation. 

1.2  Previous  Related  Investigations 

A  closed  volume  liquid-metal  combustor  utilizing  Li  and  SF6  as  fuel  and  oxidant  has  been 
identified  as  an  energy  source  for  undersea  propulsion  applications.  The  liquid-metal  combustion 
system  is  also  known  as  the  SCEPS  system  (Stored  Chemical  Energy  Propulsion  Systems),  e.g. 
see  Hughes  et  al.  (1983)  and  can  be  operated  independent  of  environmental  conditions.  The 
reaction  of  Li  and  SF6  follows  the  stoichiometry: 

8  Li  +  SFg  — »  6  LiF  +  L^S  (1) 

The  combustion  products  (LiF  and  L^S)  are  soluble  in  liquid  lithium,  yielding  two  immiscible 
liquids:  fuel-rich  (or  metal  rich)  and  product-rich  (or  salt  rich)  liquids.  The  density  of  the  product- 
rich  liquid  is  heavier  than  the  fuel  rich  liquid  (and  lithium  as  well),  allowing  for  constant  volume 
operations.  A  closed  system  operation  based  on  liquid-metal  combustion  has  been  developed  for 
undersea  propulsion  applications  (Hughes  et  al.,  1983  and  Parnell,  1989). 

As  discussed  earlier,  two  system  configurations  were  proposed,  e.g.  see  Faeth,  et  al.  (1978) 
and  Hughes,  et  al.  (1983),  the  submerged  jet  and  wick-type  or  reactive  heat  pipe  combustors.  A 
number  of  investigations  relevant  to  the  submerged-jet  operation  had  appeared  in  the  literature;  for 
example,  see  Loth  and  Faeth  (1989)  and  Parnell  (1989)  and  references  cited  therein.  Much  less 
work,  however,  was  reported  for  the  wick-type  liquid-metal  combustion,  despite  the  advantages 
over  the  submerged  jet  as  an  energy  source  for  propulsion  applications.  The  wick-type  liquid 
metal  combustion  yields  a  gaseous  diffusion  flame  near  the  wick  surface  as  identified  by  Blakeslee 
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(1977).  The  reaction  products,  however,  condensed  as  they  convected  and  diffused  to  a  low 
temperature  region  and  eventually  dissolved  in  the  liquid  bath  to  form  two  immiscible  liquids  (fuel- 
rich  and  product-rich  liquids).  A  reactive  heat  pipe  concept  based  on  the  wick  operation  was  later 
developed  and  reported  by  You,  et  al.  (1977),  Alstadt  and  Faeth  (1977),  You  and  Faeth  (1977). 
Faeth,  et  al.  (1978)  and  Groff  and  Faeth  (1978b).  The  fundamentals  of  wick-type  liquid-metal 
combustion  are  not  well  understood;  for  example,  the  presence  of  non-condensible  gases  was 
found  to  drastically  reduced  the  fuel  mass  burning  rate  (Blakeslee,  1977).  The  mechanism, 
however,  was  not  understood.  A  better  understanding  of  the  mechanism  can  provide  design 
guidance  to  reduce  the  non-condensible  gas  effects  and  provide  guidance  in  system  preparation. 

1.3  Specific  Objectives 

As  discussed,  the  wick-type  liquid  metal  combustion  has  intrinsic  merits  as  an  energy  source 
for  surface  and  under-sea  propulsion  systems  and  the  wick  configuration  can  provide  an  optical 
access  to  study  the  combustion  of  Li  and  SF6-  The  specific  objectives  of  this  project  are:  (1)  to 
document  the  flame  structure  of  wick-type  Li  and  SFg  combustion,  (2)  to  obtain  a  better 
understanding  of  wick-type  liquid-metal  combustion,  for  example  on  the  effects  of  fuel 
compositions  on  mass  burning  rates  and  (3)  to  develop  prediction  methods  for  wick  combustion 
and  to  assess  model  predictability. 
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II.  THEORETICAL  CONSIDERATION 

The  basic  phenomena  of  liquid-metal  combustion  were  discussed  in  earlier  studies;  for 
example,  the  flowfield  of  submerged  gaseous  oxidant  jets  in  liquid  metal  was  discussed  by  Avery 
and  Faeth  (1975),  Chen  and  Faeth  (1983)  and  Loth  and  Faeth  (1989)  and  the  wick-type  liquid- 
metal  combustion  was  studied  by  Blakeslee  (1977)  and  Faeth,  et  al.  (1978).  The  combustion  of  a 
wick  flame  was  found  essentially  a  gaseous  diffusion  flame  (Blakeslee,  1977)  in  which  the  fuel 
consumption  rate  was  limited  by  the  transport  in  the  gas  phase.  To  model  the  wick  combustion  of 
Li  and  SFg,  thermochemical  properties  of  reactants  and  products  were  identified  and  the  transport 
process  was  modeled. 

2. 1  Thermodynamic  Properties  of  L^S 

Thermochemical  properties  of  reactants  and  products  involved  in  Li-SFg  combustion  are 
generally  available  from  JANAF  Tables  (Chase  et  al.,  1985),  except  for  L^S  which  were  taken 
from  Groff  (1976).  This  was  done  by  curve  fitting  L^S  data  reported  by  Groff  (1976)  to  the 
format  employed  by  the  NASA  code.  The  gas  phase  property  was  estimated  from  statistical 
thermodynamics  based  on  spectroscopic  data.  Li2S  was  assumed  a  simple  harmonic  oscillator  and 
the  frequency  was  estimated  from  similar  molecules;  in  specific,  H2N,  H2O,  H2S,  and  Li20.  The 
liquid  and  solid  state  properties  were  estimated  from  the  properties  of  Na2S,  Li20,  and  Na20. 
Groff  (1976)  should  be  cons  Ited  in  estimating  L^S  properties.  To  utilize  the  NASA  equilibrium 
program  (CEC),  we  need  to  provide  specific  heat,  enthalpy  and  entropy  of  L^S  for  inclusion  in  the 
data  bank.  This  was  done  by  curve  fitting  the  L^S  data  of  Groff  (1976)  to  the  format  employed  in 
the  NASA  code. 

The  specific  heat  correlations  used  in  the  NASA  Code  (CEC  Program;  Gordon  and  McBride, 
i976)  assumed  the  following  form: 


Q 

^  =  ai  +  a2  T  +  a3  T2  +  a4  T3  +  35  T4 


(2) 
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The  coefficients  ai,  a2,  &3,  &4  and  as  were  obtained  by  curve-fitting  the  specific  heat  data  reported 
by  Groff  (1976).  A  least  squares  fit  of  fifteen  temperatures,  i.e.,  300  K  to  1600  K  at  100  K 
increments  and  at  the  L^S  normal  melting  point  (1645  K),  was  made  for  solid  phase  specific 
heats,  and  25  temperatures  (1645  K  and  1700  to  4000  K  at  100  K  increment)  and  41  temperatures 
(1000  to  5000  K  at  100  K  increment)  were  used  for  the  liquid-phase  and  gas-phase  specific  heats, 
respectively.  The  obtained  correlation  coefficients  are  summarized  in  Table  1. 


Table  1.  Coefficient  for  U2S  Specific  Heat  Correlations 


Coefficients 

Sfliid 

Liquid 

Gas 

ai 

0.56732404E+01 

0.969661 47E+01 

0.70533838E+01 

32 

0.11686870E-01 

0.11 4094 15E-02 

0.50293408E-03 

33 

-0.19245575E-04 

-0.5669227  IE-06 

-0.2277 1507E-06 

34 

0.16212627E-07 

0.12559047E-09 

0.466488 19E- 10 

35 

-0.476547 15E- 11 

-0.10217635E-13 

-0.35727596E-14 

36 

-0.61355396E+05 

-0.57484191E+05 

-0.23660 173E+05 

37 

-0.2640986 1E+02 

-0.44900746E+02 

-0.11219699E+02 

Two  additional  coefficients,  a6  and  a7,  are  also  summarized  in  Table  1  which  were  used  in  the 
enthalpy  and  entropy  correlations: 

^=ai+  ^T+  yT2+  ^T3  +  ^T4*  ^  (3) 

|  =  ailnT  +  a2T+  ^T2-^  ^T3+  ^T4^  a7  (4) 

The  two  coefficients,  a$  and  a7,  were  determined  from  the  enthalpy  integration  (a$)  and  from 
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Clausius-Claperyon  equation  (aj).  In  arriving  at  the  above  values,  the  gas-phase  enthalpy  of 
formation  was  determined,  from  statistical  mechanics,  -178.06  kJ/mole  at  298.15  K  and  0. 1  MPa. 
The  liquid-phase  enthalpy  of  formation  was  accordingly  determined  -451.55  kJ/mole  at  298.15  K. 
The  solid-phase  enthalpy  of  formation  was  determined  -492.91  kJ/mole  at  298.15  K,  based  on 
heat  of  fusion  49.4  kJ/mole  and  liquid-phase  enthalpy  -337.78  kJ/mole  at  the  normal  melting  point. 


2.2  Flow  Analysis 

The  wick  flame  was  assumed  a  diffusion  flame  with  infinitely  fast  reactions  and  boundary- 
layer  approximations  were  invoked.  The  major  assumptions  are 

•  Steady-state  two-dimensional  boundary  layer  flow. 

•  Chapman  gas,  unity  Lewis  number,  constant  Prandtl  number  and  binary  diffusion  system. 

•  No  radiative  heat  loss  and  negligible  viscous  dissipation. 

•  Infinitely  fast  chemical  reactions  and  local  equilibrium  established  in  the  flow 

•  Uniform  interface  condition. 

Based  on  the  above  assumptions,  a  conserved  scalar  approach  can  be  applied  (Williams,  1985). 
The  governing  equations  of  conservation  of  mass,  momentum,  energy  and  species  are 


Continuity 

d(pu)  3(pv)  n 
dx  +  dy 


(5) 


Momentum 


3u  9u  d  r  9u\  .  . 

pu^+pv^  =  ^[^J  +  g(p~'p) 


Mixture  Fraction 

dZ  9Z  9  (  )Z 

pu^+pv37^(pr^ 


(6) 


(7) 
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where  the  coordinate  system  is  summarized  in  Fig.  1  and  the  mixture  faction,  Z,  is  a  Shvab- 
Zeldovich  variable,  combining  the  energy  and  species  equations.  The  energy  and  species 
equations  are 


Energy 


(8) 


(9) 


Figure  1.  Coordinate  System 
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where  h  is  the  total  enthalpy,  including  enthalpy  of  formation,  and  Y;  is  the  mass  fraction  of 
species  i  and  CDi  is  the  source/sink  term  due  to  chemical  reactions.  The  governing  equations  were 
further  cast  into  a  dimensionless  form  in  terms  of  stream  function  and  similarity  variables  (Chen 
and  Faeth,  1982).  The  transformed  equatiors  are 


Momentum 


Oq  ypoo  poo 


+  3  ff  ’  -  2  (f  )2  + 


J 


P°o/p  ~  1  _  4  £  /f .  9f' 
Poo/p  w  '1  l  dt, 


(10) 


Mixture  Fraction 

a  f  p  ji  i  az)  az  afaz 

an  [poo  poo  Pr  an  J  an  q  as  as  dn 


(id 


The  continuity  is  automatically  satisfied  in  the  stream  function  formulation.  The  momentum 
equation  employs  a  dimensionless  stream  function: 


f(5.n)  = 


¥ 

4vooC(^H)i/4 


(12) 


where 


c  _  (Poq/Pw  -  1)^  1/4 

4v2 

oo 

and  the  dimensional  stream  function  was  defined  as 


d\\f  _  pu  d\\f  _  pv 

^  "pi’  ^  =  ‘pl 


(13) 


(14) 
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The  similarity  variables  are  defined  in  the  following 

\  =  cx'1/4  fy  (— )  dy;  (15) 

n  Jo  poo 

The  boundary  conditions  were  specified  to  satisfy  the  non-slip  wall  and  quiescent  ambient 
conditions: 

tj  =  0;  f  =  fw,  f  =  0,  Z  =  1,  and 

TJ  =  oo;  f  =  0,  Z  =  0.  (16) 

To  describe  the  interface  thermodynamic  state,  a  single  equation  was  obtained  for  Z  =  1.  The 
equation  is  shown  below 


where  Le,  L  and  Cpw  are  respectively  the  Lewis  number,  enthalpy  of  gasification  and  specific  heat 
at  the  interface.  Equation  (17)  was  obtained  considering  the  conservation  of  mass  and  energy  at 
the  interface,  i.e. 


Mass 


pw  Vw  (1  -  Yf,w)  —  p  D 


(18) 
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Energy 

pwvwL  =  k(|^  (19) 

The  significance  of  Eq.  (17)  is  that  the  interface  condition  is  dictated  by  the  flame  structure  within 
the  context  of  the  conserved  scalar  approach.  Eq.  (17)  yields  an  interface  fuel  mass  fraction  is  less 
than  one  under  normal  condition.  This  can  be  seen  as  the  fuel  vapor  decreases  and  the  temperature 
increases  as  Z  is  decreased  (or  moving  away  from  the  interface).  It  is  noted  that  the  interface  fuel 
vapor  mass  fraction  may  approach  unity  near  the  mixture  thermodynamic  critical  point  at  that 
condition  the  interface  fuel  mass  fraction  is  equal  to  one. 

2.3  Numerical  Methods 

As  stated  earlier,  local  equilibrium  was  assumed  for  the  flow.  The  state  relationship  was 
constructed  from  the  equilibrium  calculation.  The  equilibrium  calculation  employed  the  NASA 
equilibrium  code,  i.e.  NASA  CEC  Program  by  Gordon  and  McBride  0976).  The  CEC  program 
utilized  a  Newton-Raphsors  iterative  scheme  in  minimizing  the  Gibbs  free  energy  to  obtain 
equilibrium  species  and  temperature. 

The  flow  calculation  employed  a  finite  difference  scheme  (Keller  box  method)  to  obtain 
numerical  solutions  to  the  governing  equations,  Eqs.  10  and  11,  at  specified  boundary  and 
interface  conditions,  Eqs.  16  and  17.  The  numerical  code  was  modified,  allowing  for  wall 
blowing  effects,  from  the  computer  program  developed  by  Chen  and  Faeth  (1982).  The  modified 
program  was  tested  satisfactory  for  non-reactive  flows  (i.e.,  evaporative  cooling)  before  adopting 
it  for  reactive  flow  computations. 
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III.  EXPERIMENTAL  CONSIDERATION 
3.1  Apparatus 

It  is  known  that  Li-SF6  liquid-metal  combustors  were  operated  at  sub-atmospheric  pressures. 
Experiments  were  conducted  in  the  vacuum  chamber  shown  in  Fig.  2  which  was  modified  from  a 
commercial  product  Vacuum/Atmospheres  Model  HE- 133-5.  The  chamber  was  constructed  with 
corrosion-resistant  aluminum  alloy  (Type  6061),  having  an  interior  volume  0.43  m3  (30"  I.D.  x 
36"  long)  and  equipped  with  an  ante-chamber  (0.044  m3;  12"  I.D.  x  24"  long). 

The  chamber  interior  design  was  illustrated  in  Fig.  3  in  which  a  wick  burner  was  placed 
inside.  A  stainless  steel  wick  burner  (316  Stainless  Steel  and  100x100  mesh)  was  placed  inside 
the  chamber.  Different  burner  diameters  (32  mm  for  Li-SF6  and  25  mm  for  other  experiments)  and 
burner  heights  (25  to  89  mm)  were  used  in  the  experiments.  The  wick  screen  was  spot-welded  on 
two  stainless  steel  caps.  An  ANSI  K-type  (chromel-alumel)  sheathed  (stainless  steel  sheath) 
thermocouple  was  used  to  monitor  the  liquid  bath  temperature  and  a  thermocouple  indicator 
(Omega  DP- 80)  was  used  for  temperature  readout  The  thermocouple  output  was  also  interfaced 
with  a  Mac  IIx  work  station  using  an  A/D  converter  (National  Instrument  MIO-16-9).  A  stainless 
steel  plate  was  placed  above  the  burner  to  quench  the  hot  plume  rising  from  the  wick  burner.  A 
well  insulated  (Cole  Palmer  N-031 16-40)  brass  heater  block  (25  mm  in  diameter  and  60  mm  long) 
heated  by  nichrome  wires  (AWG  Gauge  20)  was  placed  below  the  wick  burner  and  the  heating 
was  controlled  by  a  Variac  (Fisher  Model  09-521-22)  and  monitored  by  a  Watt  meter  (Ohio 
Samitronics  EM-5).  It  required  150  W  to  heat  the  lithium  inside  the  wick  burner  to  a  molten  state 
(around  400  °C)  prior  to  combustion  experiments.  Coiled  nichrome  wires  (AWG  Gauge  20)  with 
a  small  lithium  solid  were  used  as  the  ignitor.  The  ignitor  wire  was  heated  by  a  Variac  (Research 
Lab  Vac  30-12). 


Ante-chamber 


Ante-chamber 


To  DC  Powar  Supply 
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Figure  3.  Schematic  of  the  Burner  Arrangement 

The  gas  supply  to  the  chamber  was  illustrated  in  Fig.  4.  To  prepare  the  wick  burner,  the 
chamber  was  first  evacuated  to  <0.2  torrs  using  a  chemical  pump  (Alcatel  Model  2020CP1)  and 
back  filled  with  the  zero-gas  grade  argon  gas  (Matheson  99.998%  purity)  to  atmospheric  pressure. 
A  "purge"  argon  gas  (Air  Products  zero-gas  grade,  99.99%  purity)  was  also  plumbed  to  the 
vacuum  chamber  intended  for  fast  delivery  of  argon  gas  in  case  of  emergencies.  The  pressure  of 
the  purge  gas  was  regulated  slightly  above  0.1  MPa  (1  atm).  The  chamber  pressure  was  monitored 
by  two  vacuum  gauges  (Wika,  -30730"  Hg)  and  one  pressure  transducer  (Wallace  &  Tieman 
Model  66-100,  0.005  psia  precision).  High  purity  SF6  (Matheson,  Instrument  Purity  99.99% 
purity)  was  used  as  oxidant  and  technical  grade  7Li  (99.9  %  purity.  Lithium  America;  12.7  mm 
lithium  rods  packed  in  mineral  oil  or  argon  gas)  was  used  as  fuel. 

The  lithium  sample  was  cleaned  with  hexane  and  methanol  and  carefully  prepared  to  obtain 
contamination  free  samples  for  combustion  experiments.  Prior  to  each  experiment,  the  wick 
burner  was  cleaned  with  methanol  and  baked  at  low  pressures  (inside  the  vacuum  chamber)  to 
remove  contaminants.  The  entire  chamber  (including  the  ante -chamber)  was  evacuated  to  <0.2  torr 


Pressure 

Gauge 


Pressure 

Transducer 
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Figure  4.  Gas  Supply 

and  back-filled  to  one  atmospheric  pressure  with  argon  gas.  This  step  was  repeated  at  least  once. 
The  inner  door  of  the  ante-chamber  was  then  closed  manually  (through  the  glove  port).  An 
untreated  lithium  sample  (covered  with  mineral  oil)  was  placed  inside  the  ante-chmaber  through  the 
outer  door  of  the  ante-chamber.  The  ante-chamber  was  then  evacuated  and  back-filled  to  one 
atmospheric  pressure  with  argon  gas.  The  lithium  sample  was  moved  to  the  main  chamber  from 
the  ante-chamber  through  the  inner  door  of  the  ante-chamber.  Reagent  grade  Hexane  (99.9  % 
purity)  was  used  to  remove  mineral  oil  from  the  lithium  sample  and  methanol  (Reagent  Grade;  99.9 
%  purity)  was  used  to  remove  hexane  and  oxidized  layer  from  the  lithium  sample.  Following  the 
preparation  procedure  outlined  above,  cleaned  lithium  samples  were  obtained.  The  cleaned  lithium 
sample  was  placed  inside  the  wick  burner.  The  chamber  was  evacuated  and  the  lithium  sample 
was  heated.  The  lithium  sample  was  heated  to  400  °C  (or  above)  to  wet  the  wick  surface.  The 
chamber  was  then  back-filled  with  SFg  for  subsequent  combustion  experiments.  The  combustion 
process  was  monitored  and  recorded  using  a  VHS  format  video  camera  (Panasonic  WV31 10-DX) 
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for  subsequent  review  and  analysis.  A  35  mm  camera  (Canon  AE-1)  was  used  to  take 
photographs  of  the  video  directly  off  a  TV  screen. 

The  optical  arrangement  of  the  diagnostics  is  shown  in  Fig.  5  in  which  a  top  view  of  vacuum 


Compartment  Viewing  Port 


Figure  5.  Experimental  Set  Up  for  Wick  Combustion  Diagnostics 
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chamber  was  illustrated.  A  three  dimensional  traversing  mechanism  which  consists  of  linear 
bearings  and  uni-slides  was  designed,  acquired  (from  commercial  products)  and  assembled  with 
in-house  made  components.  The  design  provides  a  fixed  optical  axis  (focusing  lenses  and 
spectrometer)  for  laser  induced  fluorescence  (LIF)  measurements;  the  variation  in  vertical  heights 
was  accounted  for  by  a  pair  of  periscopes.  Similar  to  that  employed  in  the  layout  of  LIF 
measurements,  another  traversing  mechanism  was  also  designed  and  assembled  for  laser  Doppler 
velocity  measurements  using  a  back  scattering  mode  operation. 

3 . 2  LIF  Thermometry 

To  develop  the  laser  induced  fluorescence  (LIF)  thermometry  using  7Li2  (Lithium  America, 
99.9  %  purity)  as  fluorescence  species,  experiments  were  conducted  with  the  apparatus  shown  in 
Fig.  6.  The  major  components  of  the  apparatus  include  a  heat-pipe  oven,  an  argon-ion  laser, 
optics  and  photodetectors.  The  heat  pipe  was  configured  with  five  arms.  Each  arm  has  an  optical 
window  at  its  end.  The  heat  pipe  was  operated  to  maintain  a  uniform  inside  temperature  and  was 
used  for  LIF  thermometry  calibration  experiments.  The  oven  temperature  was  determined  from 
thermocouples  (ANSI  K-type)  located  at  the  center  of  the  heat  pipe.  In  the  experiments,  the  central 
portion  of  the  heat  pipe  was  electrically  heated  and  the  far  end  was  water  cooled  to  protect  the 
windows  from  attack  by  Li  vapor.  In  preparing  calibration  experiments,  lithium  was  placed  inside 
the  heat  pipe  and  followed  by  evacuation  and  heating  to  remove  contaminants  from  lithium.  Argon 
gas  was  introduced  to  confine  Li  vapor  in  the  central  portion  of  the  heat  pipe,  to  protect  the 
window  from  direct  contact  with  lithium  vapor  and  to  control  the  system  pressure.  Adjusting 
heater  and  vacuum  pump  operation,  steady-state  conditions  of  desired  temperature  and  pressure 
can  be  reached  for  LIF  thermometry  experiments. 

An  argon  ion  laser  (Coherent  CR-6,  6  W)  was  used  to  excite  7Li2-  The  coincidence  lines  of 
argon  ion  laser  output  and  7Li2  spectrum  are  well-documented  (e.g.,  Luh,  Bahns  and  Stwalley, 
1988),  the  values  are  summarized  in  Table  2.  In  the  present  study,  the  single  line  excitation  was 


Table  2.  Fluorescence  Lines  of  7Li2  Molecule  Excited 
by  Argon  Ion  Laser  (B-X  Transitions) 


Laser  Wavelength 

B-X 

LIF  Wavelength 

(nm) 

fv'J’  -  v"J") 

fnm) 

457.9 

(9,9  -2,8) 

(9,5  -  2,5) 

533.9-567.2 

550.3-564.9 

472.7 

(8,32-3,31) 

550.0-568.6 

476.5 

(4,24-1,24) 

(9,61-2,62) 

(11,34-5,34) 

525.1-551.0 

582.7-586.4 

551.0-588.9 

488.0 

(2,31-1,30) 

(6,45-3,46) 

(4,73-0,74) 

(9,19-6,18) 

(12,36-7,36) 

514.7- 540.7 
548.4-559.6 
552.3-564.7 

566.8- 568.1 
584.1-591.3 

514.5 

(8,61-7,61) 

(3,30-5,29) 

(9,38-9,37) 

574.6-582.1 

Halt  Pip*  Oven 


Figure  6.  LIF  Thermometry  Calibration  Experimental  Set  Up 
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employed  In  the  discussion  herein  the  temperature  was  obtained  from  a  single  line  excitation  (488 
nm),  assuming  that  the  two  closely-spaced  fluorescence  lines  excited  by  a  single-laser  line  have 
constant  collection  volume,  detection  efficiency  and  same  quenching  effects.  The  fluorescence 
signal  was  collected  at  a  right  angle  from  the  excitation  laser  beam.  A  0.85  m  double  spectrometer 
(SPEX-1404  )  was  used  to  obtain  the  fluorescence  spectrum.  The  height  of  the  spectrometer 
entrance  slit  to  was  20  mm  and  the  width  was  adjusted  in  the  range  20  to  50  p.m  to  obtained  good 
spectral  resolutions.  The  measured  fluorescence  spectrum  was  analyzed,  comparing  with  the 
theoretical  spectrum.  A  Xe-lamp  was  used  as  reference  to  ensure  the  wavelength  assignments. 

To  examine  the  fluorescence  spectrum  in  a  quantitative  manner,  the  theoretical  frequency  and 
corresponding  intensity  distribution  must  be  calculated  from  the  theory  to  confirm  the  assigned 
fluorescence  spectrum.  The  molecular  constants  of  B^Iu  and  X1!^  states  of  7Li2  were  taken  from 
Hessel  and  Vidal  (1978)  to  calculate  theoretical  fluorescence  line  frequencies  from  specified  excited 
states: 

v  =  T(v',J')  -  T"(v",J")  +  AT(v',J')  (20) 


where  v  is  the  fluorescence  frequency  resulting  from  the  difference  in  term  energies  between  the 
excited  state  T(v’,J')  and  the  ground  state  T"(v",J").  The  v  and  J  are  the  vibrational  and  rotational 
quantum  number  of  each  respective  state.  The  term  energy  of  the  excited  B  level  (v',J')  is 

T(v  ’,  J’)  =  I  A’j:  (v'+  t  y  [  J’  (J’+l)  -  1  ]  J  (21) 

ij  J 

and  the  corresponding  term  value  of  the  ground  state  is  given  by 
T(v",  J")  =  £A"jj  (v"4)>[J"  (J"+l)]j 

y 


(22) 


19 


The  lambda-doubling  AT(v',J')  of  Btflu  state  is  as  follows 

AT(v’J')  =  8  (  I  a'lk  (v'+  \  )i  (J*  (J'+l)  -  1)  j  )  (23) 

lk  x 

where  8  =  0  for  the  Q  series  (AJ  =  0)  and  8  =  1  for  P  and  R  series  (AJ  =  ±  1).  The  Dunham- type 
molecular  constants  are  denoted  as  A'jj  A"jj  and  a'^  (Hessel  and  Vidal,  1978).  The  theoretical 
intensities  of  fluorescence  spectrum  were  calculated  following  Luh  (1985)  in  which  a  computer 
program  known  as  Intensity  (Zemke  and  Stwalley,  1984)  was  used  to  sofre  the  Schrodinger 
equation  for  theoretical  transition  frequencies  and  the  corresponding  intensities. 

In  Li2  LEF  thermometry  calibration,  the  488.0  nm  output  of  the  Argon-ion  laser  was  used  as 
the  excitation  source.  The  resultant  fluorescence  spectrum  is  a  collection  of  7Li2,  6Li2  and  7Li6Li 
molecules.  The  fluorescence  of  7Li2  resulted  in  two  stronger  series  from  the  excited  states  B(2,3l) 
and  B(6,45).  Die  ratio  of  spectrally  resolved  fluorescence  lines  from  B(2,31)  and  B(6,45)  was 
used  to  obtain  the  temperature  information.  Two  models,  linear  fluorescence  and  saturated 
fluorescence,  were  considered  in  analyzing  the  spectroscopy  data.  The  measured  fluorescence 
signal  can  be  expressed  as 

0  =  PQcVc(^)NjIL  (24) 

dS2 

where  |3  denotes  the  collecting  efficiency,  Qc  die  solid  angle,  Vc  the  collection  volume,  (•— )  the 

dQ 

scattering  cross  section,  Nj  the  number  density  at  the  excited  state  and  II  the  laser  flux. 

Linear  Fluorescence.  A  three-level  model  was  employed  to  describe  the  fluorescence  signal 
resulting  from  transitions  due  to  laser  excitation.  The  linear  fluorescence  is 


O-PQcVch  vjk  Ajk  Nj 


(25) 
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where  h  is  the  Planck's  constant,  Vjk  and  Aj*  are  the  frequency  and  Einstein  coefficient  of  the 
florescence  transition  between  the  excited  state  j  and  ground  state  k.  The  population  of  the  excited 
state  under  the  steady  state  assumption  can  be  expressed  as  a  function  of  absorption  Wjj,  Einstein 
coefficient  A  and  quenching  rate  Q: 


Ni"Wji 


W; 


!L 


+  A  +  Q 


(26) 


The  absorption  Wjj  is  related  to  Einstein  coefficient  Ay: 


Bjj  iL(vii) 
Jt  c  Av 


(27) 


where  Av  is  the  spectral  width  of  the  laser  line  and  the  absorption  coefficient  By  is  defined  as 


By- 


8;ch  vy 


3  AiJ- 


(28) 


The  above  equation  was  obtained  assuming  negligible  quenching  effects  (Q  =  0)  which  is  typical  at 
low  pressures  and  negligible  stimulated  emission  (Wjj  =  0)  compared  to  the  absorption  Wy.  In 
fact,  the  stimulated  emission  can  be  eliminated  when  the  florescence  signal  was  collected  at  a  right 
angle  from  the  laser  beam.  It  was  further  assumed  that  the  population  follows  a  Boltzmann 
distribution,  yielding  a  linear  equation  for  fluorescence 


„  =  (ivMAjiM!)  Vjk  {Ak 

16kJc  Av  J  ZAi’tr 


x  KT  gi_  ,hcEi 


(29) 


Where  Av  is  the  laser  spectral  width,  Nj  the  total  number  density  of  laser  coupled  states,  gj  the 
degeneracy  of  the  excited  state  (2Jj  +  1),  Zjot  the  total  partition  function  and  Ei  the  energy  of  the 
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initial  state. 

Considering  two  fluorescence  lines  resulting  from  two  different  series  excited  by  the  same 
laser  line,  one  can  extract  the  temperature  from  the  ratio  of  the  two  fluorescence  signal  Otin!  and 
^lin^ 


T  = - hc_AE -  (30) 

£  ln  r^linl  A2  V2  g2  Il(V2)  ^  B  2 
t^li.12  Al  V!  g,  IL(V,)  J  B1 


where  AE  is  the  ground  state  energy  difference  (cm-1)  between  the  two  fluorescences  lines  of 
different  series.  Since  the  laser  intensity  is  likely  to  vary  over  the  two  excitation  frequencies,  Vi 
and  V2,  the  difference  in  Il(Vi)  and  I1XV2)  must  be  considered  when  the  linear  fluorescence  model 
is  employed  for  temperature  measurements. 

Saturated  Fluorescence.  At  saturated  states,  the  fluorescence  signal  is  independent  of  laser 
power  and  absorption  coefficients.  The  saturated  fluorescence  signal  is  (Lucht,  1982): 


<I>sat  «  Ajk  vjk  exp 


h  c  EjN 
K  T 


(31) 


where  gj  and  gj  are  the  degeneracies  of  the  ground  and  excited  states,  respectively.  The  above 
equation  also  assumes  negligible  quenching  effects.  Taking  the  ratio  of  two  saturated  fluorescence 
line  signals,  Osatj  and  Osat<2,  one  obtains 


_ h  c  AE _ 

f  S2i  g2i  \ 

frsatl  A2  V2  g2i  +  M2  i 

3)sat2  Ai  Vi  &li  glj 

V  gli  +  glj 


02) 


where  temperature  is  expressed  in  terms  of  ground  state  energy  difference  (AE),  saturated 
fluorescence  intensities  (Osat),  frequencies  (v),  Einstein  coefficients  (A)  and  degeneracies  (g). 
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IV.  RESULTS  AND  DISCUSSION 
4. 1  Experimental  Results 

LIF  Thermometry.  The  fluorescence  spectrum  of  7Li2  excited  by  an  Argon  ion  laser  at  488 
nm  was  analyzed.  The  assignments  of  the  spectrum  were  confirmed  with  theoretical  calculations  in 
terms  of  fluorescence-line  wavenumbers  and  relative  intensities.  The  spectral  assignments  are 
summarized  in  Tables  3  and  4  where  X  is  the  wavelength  in  A  and  CD  is  the  vacuum  wavenumber. 
The  comparison  of  relative  intensity  between  experimental  data  and  theoretical  predictions  was 
shown  in  Figs.  7  and  8.  The  uncertainty  of  intensity  measurements  can  be  improved  by  reducing 
the  spectrometer  scanning  speed  and  spectrometer  entrance  slit  width.  Setting  the  scanning  speed 
at  5  A/s  and  the  entrance  slit  at  25  (im,  the  measured  fluorescence  intensity  at  507  "C  was  in  good 
agreement  with  theoretical  prediction,  cf.  Fig.  7.  The  measured  fluorescence  ratios  of  three 
different  transition  series  in  the  temperature  range  850  to  1050  K  were  compared  to  theoretical 
calculations.  The  experimental  data  were  bounded  by  the  calculations  employing  linear  and 
saturated  fluorescence  models,  Figs.  9-11.  As  can  be  seen  the  comparison  is  encouraging,  the 
experimental  data  were  found  bounded  by  linear  fluorescence  and  saturated  fluorescence  model 
calculations.  Refinements,  however,  are  needed  to  implement  LIF  thermometry  techniques  to 
combustion  studies.  In  specific,  the  effects  of  excitation  laser  intensity  on  measured  fluorescence 
signals  need  to  be  evaluated. 

Wick  Combustion.  Prior  to  Li-SFg  combustion  experiments,  the  vacuum  chamber  was 
tested  with  combustion  experiments  of  ethanol  and  hexane  burning  in  the  air.  A  vertical  cylindrical 
wick  (25  mm  I.  D.  and  50  to  89  mm  long)  in  which  nonflammable  fiberglass  was  packed  was 
placed  inside  the  chamber.  The  fiberglass  was  used  to  retain  the  liquid  fuel  inside  the  wick. 
Reagent  grade  ethanol  (99.9  %  purity)  or  hexane  (99.9  %  purity)  was  used.  The  wick  combustion 
yields  a  non-luminous  gaseous  diffusion  flame  for  ethanol  at  low  pressures  (0.03  and  0.05  MPa) 
and  a  luminous  flame  for  ethanol  and  hexane  at  atmospheric  pressure  (0.1  MPa),  e.g.  see  Figs.  12 
and  13  for  ethanol  at  0.03  and  0.1  MPa  and  Fig.  14  for  hexane  at  0.1  MPa.  The  combustion 
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process  can  be  divided  in  three  stages:  ignition,  stable  burning  and  extinction.  After  the  ignitor 
was  activated,  a  stable  axisymmetric  diffusion  flame  around  the  wick  was  established  following  a 
transient  ignition  stage.  During  the  ignition  stage,  the  flame  spread  over  the  wick  surface,  it 
usually  started  from  the  upper  portion  of  the  wick  and  propagated  down  to  the  leading  edge  of  the 
cylinder.  Near  the  fuel  burn-out,  the  flame  became  unstable.  Local  extinction  from  the  leading 
edge  of  the  wick  was  seen,  yielding  an  oscillating  flame  surface  in  longitudinal  direction  until 
extinction  was  observed. 


Table  3.  Assignment  of  7Li2  B(2,31)  Fluorescence  Excited  at  488  nm 
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Table  4.  Assignment  of  7Li2  B(6,45)  Fluorescence  Excited  at  488  nm 
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Relative  Intensity  of  7Li2  B(2/31)-X(v",30) 


Figure  7.  Relative  LIF  Intensity  of  7Li2  B(2,31)  Excited  at  488  nm. 
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Figure  9.  LIF  Thermometry  Calibration  of  7Li2;  B  (2,  31)-X(0,  30/32)/B(6,45)-X(2, 44/46) 
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Figure  10.  LIF  Thermometry  Calibration  of  7Li2;  B  (2,  31)-X(5,  30/32)/B(6,45)-X(8, 44/46) 
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Thermometry  Calibration  of  7Li2;  B  (2,  31)-X(6,  30/32)/B(6,45)-X(8,44/46) 
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Spatial  and  temporal  periodicities  were  observed  in  wick  diffusion  flames.  The  oscillation 
frequency  was  estimated  in  the  range  5  to  8  Hz.  The  oscillation  frequency  of  ethanol  diffusion 
flames  decreased  when  the  system  pressure  was  increased.  Ethanol  and  hexane  were  found  to 
have  similar  oscillation  frequencies  at  0.1  MPa.  The  spatial  wavelength  (e.g.  30  to  80  mm  in  the 
present  experiments)  increased  as  the  system  pressure  was  decreased.  It  is  also  noted  that  the 
oscillation  in  wick  diffusion  flames  is  similar  to  flame  flickers  in  buoyant  jet  diffusion  flames 
(Chen  et  al,  1989b),  suggesting  the  existence  of  buoyancy  driven  vortical  structures  in  a  vertical 
boundary-layer  diffusion  flame.  The  flame  stand-off  distance  increased  as  the  system  pressure 
was  decreased  (around  5  mm  at  0.1  MPa  and  10  mm  at  0.03  MPa  at  a  location  25  mm  from  the 
leading  edge),  in  qualitative  agreement  with  theoretical  predictions. 

Li-SF*  Combustion.  Combustion  of  lithium  solids  (in  flake  form,  i.e.  2  to  3  mm  thickness 
and  12  mm  in  diameter)  in  SF6  (0.01  MPa)  was  successfully  accomplished  using  Nichrome  wires 
(AWG  Gauge  20)  as  an  ignitor.  Following  a  transient  stage  (ejection  of  lithium)  after  the  ignitor 
was  activated,  a  very  bright  flame  covered  the  lithium  solid.  A  lithium  "droplet"  (molten  lithium) 
war-  subsequently  formed  on  the  stainless  steel  plate  inside  the  combustion  chamber.  Lithium  was 
seen  to  wet  the  stainless  steel  (spreading  on  the  surface  in  a  dome  or  bell  shape).  The  burning  of 
the  Li  drop  continued  for  a  few  seconds  before  before  the  extinction  and  solidification  of  the  melt 
occurred.  The  sequential  events  of  Li-solid  combustion  were  summarized  in  Fig.  15. 

Wick  combustion  of  Li  and  SF6  was  studied.  In  the  experiments  conducted  up-to-date,  the 
chamber  was  filled  with  SF6  at  0.01  MPa.  The  diameter  of  the  wick  was  32  mm  and  the  height, 
excluding  the  end  caps,  was  around  25  mm.  Only  one  layer  stainless  steel  screen  was  used  to 
construct  the  wick  cylinder.  Heating  the  lithium  to  a  temperature  above  400  'C  was  necessary  to 
wet  the  wick  surface.  A  long  nichrome  wire  coil  (about  200  mm  long)  sided  with  the  wick 
cylinder  was  used  as  the  ignitor  in  the  experiment  shown  in  Fig.  16.  Photographs  shown  in  Fig. 
16  were  taken  from  the  video  at  six  different  times  after  the  ignitor  was  actuated,  i.e.  2,  4,  17,  29, 
32  and  33  seconds  (cf.  Fig.  16,  from  left  to  right  and  top  to  bottom).  Once  successful  ignition  was 
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Figure  15.  Li  Solid-SF6  Combustion  at  0.01  MPa 


achieved,  the  flame  sustained  itself  until  the  wick  was  burnt  out.  During  combustion  the  chamber 
interior  showed  a  glowing  pinkish  color.  A  bright  pinkish  flame  around  the  wick  was  recorded  in 
the  video  tape.  The  aperture  of  the  video  camera  was  partially  closed.  Immediately  after  ignition, 
the  flame  spread  over  the  wick  surface  and  dark  spots  inside  the  luminous  flame  were  observed. 
The  flame  then  covered  the  side  (wick)  and  bottom  (cap)  surfaces  of  the  cylinder,  showing  a 
"uniform"  flame  (bright)  surface.  The  flame  surface  (bright  luminous  layer)  appeared  quite  close 
to  the  wick  surface,  the  stand-off  distance  was  only  a  few  mm  away  from  the  wick  surface.  Near 
bum  out,  local  extinction  first  appeared  at  the  bottom  surface  and  the  upper  portion  of  the  wick. 
The  extinction  of  the  flame  at  the  upper  wick  surface  probably  was  a  result  of  non-wetted  wick 
surface  near  "complete  consumption"  of  lithium  fuel.  Shortly  after  local  extinction  at  the  upper 
wick  surface,  the  wick  was  burned  out. 

4.2  Theoretical  Results 

Ethanol-Air  Wick  Diffusion  Flame.  The  conserved  scalar  approach  yields  a  single  equation 
for  the  interface  condition,  Eq.  (17),  dictated  by  the  flame  structure.  One  immediate  observation  is 
that  the  interface  fuel  vapor  mass  fraction  is  less  than  one  under  normal  condition.  This  can  be 
seen  as  the  fuel  vapor  mass  fraction  decreases  and  the  temperature  increases  as  Z  is  decreased  (or 
moving  away  from  the  interface).  It  is  noted  that  the  interface  fuel  vapor  mass  fraction  may 
approach  unity  near  the  mixture  thermodynamic  critical  point.  Two  system  pressures,  0.03  and 
0.1  MPa,  were  considered  for  ethanol-air  wick  combustion.  The  interface  temperatures  were 
determined  310  K  and  336  K  and  the  interface  fuel  vapor  mass  fractions  were  0.6289  and  0.6382 
for  0.03  and  0. 1  MPa,  respectively.  The  interface  condition  was  obtained  employing  an  iterative 
procedure  to  solve  Eq.  (17).  Up  to  ten  iterations  were  required  for  the  two  pressures  considered. 
The  state  relationships  of  ethanol-air  wick  diffusion  flames  at  0.03  MPa  and  0.1  MPa  were 
illustrated  in  Figs.  17  and  18,  respectively.  The  higher  flame  temperature  at  0.1  MPa  is  a  result  of 
the  equilibrium  shift  due  to  reduced  dissociation  at  higher  pressures. 


Mixture  Fraction 


Figure  17.  State  Relationship  (Temperature  and  Density)  of  Ethanol-Air  Wick  Diffusion  Flame 


Figure  18.  State  Relationship  (Concentration)  of  Ethanol- Air  Wick  Diffusion 
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In  the  computation  reported  herein,  a  flat-plate  100  mm  wick  was  considered.  The 
computation  showed  that  similarity  solutions  were  maintained  for  the  conditions  examined 
although  non-similar  governing  equations  were  retained  in  the  analysis.  The  solution  yields  similar 
profiles  for  f,  f,  f',  Z  and  Z'  at  different  streamwise  locations  when  similarity  variables  were 
employed.  The  calculated  mass  burning  rates  yields  a  £-0.25  dependence,  typical  of  laminar 
boundary  layer  diffusion  flames.  The  calculated  mass  burning  rates  are  illustrated  in  Fig.  19.  The 
predicted  mass  burning  rates  also  showed  a  pressure  dependence  of  p.0-64 

To  illustrate  the  pressure  effects  on  the  flame  stand-off  distance  from  the  wick  surface,  the 
predicted  temperature  and  streamwise  velocity  profiles  at  the  trailing  edge  of  the  100  mm  wick 
were  shown  in  Fig.  20.  As  shown  in  the  figure,  the  stand-off  distance  increased  as  the  system 
pressure  was  decreased.  This  is  in  qualitative  agreement  with  the  experiments  conducted,  although 
the  analysis  over-estimated  the  stand-off  distance  by  50  %.  The  maximum  streamwise  velocities 
are  comparable  for  the  two  pressures  considered,  with  a  higher  maximum  velocity  at  0. 1  MPa. 
This  is  not  surprising  as  the  momentum  source  term  in  dimensionless  form  is  similar  for  the  two 
pressures  considered. 

Li-SF6  Wick  Diffusion  Flame.  Equilibrium  computation  was  made  to  estimate  the  flame 
structure  of  the  wick  combustion,  to  estimate  the  pressure  effects  on  wick  combustion  and  to 
estimate  the  fuel  composition  effects  on  lithium  combustion.  The  state  relationships  of  Li-SF6  wick 
flame  at  0.01  and  0.1  MPa,  assuming  Yy  =  1  at  the  interface,  are  illustrated  in  Figs.  21-23.  At  the 
interface  (Z  =  1),  the  fuel  vapor  was  set  at  1339  K  and  1609  K  (saturation  temperatures)  for  0.01 
and  0.1  MPa,  respectively.  The  equilibrium  computation  yields  similar  state  relationships  for  0.01 
and  0.1  MPa.  The  0.1  MPa  case,  however,  has  a  higher  temperature  (4360  K  versus  3884  K  at 
0.01  MPa)  and  a  higher  density  as  shown  in  Fig.  20.  Fuel  vapor  was  the  major  species  at  Z  =  1 
either  as  lithium  atoms  (Li)  or  lithium  dimers  (Li2),  e.g.  see  Figs.  22  and  23.  The  oxidant 
appeared  as  SF6  or  SF4  at  Z  =  0;  however,  SFg  decreased  but  SF4  increased  as  Z  increased  from 
zero.  From  the  equilibrium  prediction,  LiF  is  the  major  combustion  product.  Gas-phase  LiF 


T 


Ethanol-Air 

S  0.1  MPa 

0.03  MPa  ' 

Ethanol-Air 

'' -  0.1  MPa 

. <. 

0.03MPa^^^''  " 

Figure  19.  Predicted  vw  and  m"  of  Ethanol-Air  Wick  Diffusion  Flames 
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Figure  20.  Predicted  Flame  Structure  of  Ethanol-Air  Wick  Diffusion  Flames 
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Figure  21.  State  Relationship  of  Density  and  Temperature  of  Li-SF6  Wick 


Combustion 


Figure  23.  State  Relationship  of  Species  Concentration  of  Li-SF<j  Wick  Combustion  at  0. 1  MPa 
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appeared  over  much  of  the  mixture  fraction  space  while  liquid-  and  solid-phase  LiF  appeared  in  the 
oxidant  rich  regime.  Only  gaseous  L^S  in  the  fuel  rich  regime  was  predicted  from  the  equilibrium 
calculation.  The  S  species  appeared  as  elemental  sulfur  (S)  or  SFX  (x  =  1,2,  ..,  or  6)  in  the 
oxidant-rich  regime. 

The  flowfield  solution  of  Li-SFg  wick  combustion  at  0.01  and  0.1  MPa  was  summarized  in 
Figs.  24  and  25,  employing  the  state  relationships  shown  in  Fig.  21.  Similar  to  ethanol-air  wick 
combustion,  the  prediction  showed  that  similarity  profiles  prevailed  in  the  flow.  For  example, 
nearly  identical  mixture  fraction  profiles  (Z)  and  dimensionless  velocity  (f )  were  found  at  different 
streamwise  locations,  cf.  Figs.  24  and  25.  The  pressure  was  also  found  to  have  negligible  effects 
on  predicted  mixture  fraction  and  dimensionless  velocity  profiles  as  one  compared  the  results  of 
0.01  MPa  and  0.1  MPa.  The  peak  dimensionless  velocity  appeared  at  q  =  1,  corresponding  to  the 
maximum  temperature  location  at  which  Z  =  0.3.  The  lower  pressure  condition,  i.e.  0.01  MPa, 
has  a  higher  maximum  dimensionless  velocity  (f  =  0.40)  than  the  0.1  MPa  condition  (f  =  0.38). 
The  streamwise  flow  development  was  summarized  in  Figs.  26  and  27.  The  maximum  streamwise 
velocity  increased  from  3.8  m/s  at  ^  =  0.22  (x  =  0.022  m)  to  8.1  m/s  at  ^  =  1  (x  =  0.1  m)  for  the 
condition  at  0.01  MPa  while  the  maximum  velocity  increased  from  4.0  m/s  (£  =  0.22)  to  8.5  m/s 
=  1.0)  for  the  condition  at  0.1  MPa.;  it  is  note  that  the  0.1  MPa  condition  had  a  lower  maximum 
dimensionless  velocity  (cf.  Fig.  25). 

The  boundary  layer  thickness  grows  in  the  streamwise  direction,  roughly  following  a  x0-25 
dependence  similar  to  that  expected  for  natural  convection  laminar  boundary  layer  flows.  The 
flame  stand-off  distance  increased  as  the  system  pressure  was  decreased;  for  example,  the  stand¬ 
off  distance  at  0.01  MPa  was  predicted  three  to  four  times  the  value  at  0.1  MPa,  cf.  Figs.  26  and 
27.  The  stand-off  distance,  however,  remained  constant  over  the  entire  plate  in  the  transformed 
domain.  This  confirms  the  proper  similarity  transformation  used  in  the  analysis.  The  local  wall 
blowing  velocity  due  to  gasification  of  the  liquid  was  shown  in  Fig.  28  and  the  local  fuel  mass 
burning  rate  per  unit  surface  area  was  summarized  in  Fig.  29.  The  wall  blowing  velocity  (vw)  was 
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Figure  27.  Temperature  and  Velocity  Profiles  of  Li-SF6  Wick  Combustion  at  ^  =  1 
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determined  from  Eq.  (18)  and  the  fuel  mass  burning  rate  (m”)  was  obtained  by  multiplying  vw  with 
local  gas  density.  Both  vw  and  m"  revealed  a  x-°-25  dependence,  similar  to  the  results  obtained  for 
ethanol-air  wick  diffusion  flames,  in  agreement  with  the  classical  results  of  heat  and  mass  transfer 
in  natural  convection  laminar  boundary  layer  flows. 

Also  shown  in  Figs.  28  and  29  are  the  results  based  on  the  equilibrium  properties  of  fuel-rich 
and  product-rich  liquids  as  fuel.  The  species  composition  of  the  fuel  was  taken  from  Groff  and 
Faeth  (1978a)  for  the  condition  at  0.0035  MPa  (1250  K),  e.g.  see  Fig.  30.  The  enthalpy  of 
gasification  was  estimated  assuming  an  ideal  solution  at  its  proper  composition.  The  predicted 
mass  burning  rate  decreased  as  the  system  pressure  was  decreased  when  pure  lithium  and  fuel-rich 
liquids  were  considered;  it  is,  however,  quite  surprising  to  note  that  the  product-rich  condition  at 
0.0035  MPa  yields  a  mass  burning  rate  higher  than  that  of  the  fuel-rich  condition  at  0.0035  MPa 
and  that  of  pure  lithium  at  0.01  MPa  (cf.  Fig.  29).  The  higher  mass  burning  rate  is  a  result  of  the 
increased  heat  transfer  from  the  flame  to  the  wick  surface  as  the  flame  stand-off  distance  decreases 
when  the  initial  Li  mass  fraction  was  decreased.  The  higher  mass  burning  rate  is  also  a  result  of 
reduced  enthalpy  of  gasification  per  unit  mass  of  the  liquid,  e.g.  21228  kJ/kg  for  pure  Li  1609  K, 
21950  for  pure  Li  at  1339  K,  and  21631  kJ/kg  for  pure  Li,  19572  kJ/kg  for  fuel-rich  liquid  and 
8381  kJ/kg  for  product-rich  liquid  at  1250  K.  The  higher  mass  burning  rate  does  not  imply  a 
higher  Li  consumption  rate  since  only  little  Li  was  present  in  the  product-rich  fuel.  The  results, 
however,  suggest  that  the  wick  configuration  may  have  the  potential  to  achieve  a  higher  fuel 
utilization  than  the  submerged  jet  operation. 


Figure  31.  Li-SFe  Wick  Combustion  of  Fuel-Rich  and  Product-Rich  Fuel  at  0.0035  MPa. 
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V.  SUMMARY  AND  FUTURE  WORK 

An  experimental  and  theoretical  investigation  was  conducted  to  study  wick  combustion  of  Li- 
SF6.  A  single-line  laser  induced  fluorescence  thermometry  technique  with  U2  as  fluorescence 
species  was  developed  and  calibration  experiments  showed  promising  results  with  refinements 
identified  for  further  development.  Wick  combustion  of  ethanol  and  hexane  in  air  as  well  as  Li  and 
SF6  was  conducted  in  a  vacuum  chamber  at  atmospheric  and  sub-atmospheric  pressures.  The 
ethanol  and  hexane  wick  diffusion  flames  showed  flame  oscillations,  similar  to  flame  flickers  in 
buoyant  jet  diffusion  flames.  The  oscillation  frequency  was  estimated  in  the  range  5  to  8  Hz  at  the 
pressures  examined.  The  flame  stand-off  distance  was  found  to  increase  when  the  system  pressure 
was  decreased.  The  wick  combustion  of  Li  and  SF6  resulted  in  a  bright  pinkish  flame.  The 
luminous  zone  appeared  quite  close  to  the  wick  surface  at  the  condition  examined.  Near  the 
"complete  consumption"  of  lithium,  the  wick  was  burned  out  due  to  loss  of  lithium  as  a  heat  sink. 

On  the  analysis,  a  conserved  scalar  approach  was  employed  to  model  the  wick  flame.  A 
single  equation  was  obtain  to  describe  the  interface  condition  of  wick  combustion.  Numerical 
solutions  were  obtained  for  laminar  wick  diffusion  flames.  The  prediction  yields  similarity  profiles 
for  both  ethanol-air  and  Li-SFg  wick  diffusion  flames  although  non-similar  governing  equations 
were  retained  in  the  formulation.  The  flame  stand-off  distance  was  found  to  increase  as  the  system 
pressure  was  reduced,  in  agreement  with  theoretical  predictions.  The  effects  on  mass  burning  rates 
due  to  varying  fuel  composition  were  studied  employing  the  fuel-rich  and  product-rich  liquids  at 
0.0035  MPa  in  the  analysis.  It  was  found  the  product-rich  fuel  yields  a  much  lower  flame 
temperature  but  the  mass  burning  rates  was  not  decreased.  This  is  a  result  of  the  decrease  in  the 
flame  stand-off  distance,  yielding  an  increase  in  heat  transfer  to  the  wick  surface,  as  well  as  the 
reduced  specific  enthalpy  of  gasification  when  a  product-rich  liquid  was  used.  In  fact,  the  mass 
burning  rates  were  higher  than  the  fuel-rich  fuel  at  the  same  pressure  and  even  higher  than  the  pure 
lithium  at  0.01  MPa.  The  findings  on  fuel  mass  burning  rates  are  interestir"  since  it  suggests  that 
the  wick  configuration  may  achieve  a  higher  fuel  utilization  than  the  submerged  jet. 
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Work  is  in  progress  to  document  the  spectroscopy  of  Li-SF6  wick  diffusion  flames  and  to 
include  the  curvature  effects  in  modeling  wick  combustion.  It  should  be  noted  that  the 
spectroscopy  data  of  Li  and  SF6  combustion  are  not  presently  available  in  the  literature.  The 
current  work  represents  the  first  experimental  effort  to  obtain  such  data.  The  inclusion  of  curvature 
in  the  modeling  will  clarify  its  effects  on  fuel  mass  burning  rates  and  on  the  transport  processes. 
The  wick  flames  were  found  to  exhibit  spatial  and  temporal  periodicities  in  ethanol  and  hexane 
diffusion  flames  and  time  dependent  burning  of  Li  in  SF6.  The  modeling  will  be  extended  to 
include  the  time  dependent  effects  and  the  experiments  will  be  made  for  time -dependent  Li  mass 
burning  rate  measurements.  In  quantifying  the  thermochemical  properties  of  Li  and  SF6 
combustion,  the  feasibility  of  LIF  with  LiF  as  the  fluorescence  species  will  be  explored.  This  is 
due  to  that  the  signal  to  noise  ratio  of  Li2  LIF  may  be  deteriorated  in  the  combustion  environment 
as  a  result  of  the  strong  luminosity  observed  in  the  experiments.  Thus  it  is  desired  to  use  the  major 
combustion  product,  LiF,  as  the  fluorescence  species. 

New  laser  diagnostics  will  also  be  employed  in  the  next  year's  experiments;  for  example,  a 
gated  image  intensified  (IIT  F41 1 1  Intensifier  Tube  and  Power  Supply  @  180  ns  shutter  opening) 
CCD  Camera  (Photometrics/Thomson  CSF  TH7882CDA  CCD;  576  x  384  Pixels)  will  be  used  for 
planar  imaging  of  the  wick  Li-SF6  combustion.  Due  to  the  strong  flame  luminosity,  the  frequency 
doubled  second  harmonic  output  of  a  Nd:YAG  laser  (Spectra  Physics  GCR-4-10  Nd:YAG  Laser 
and  PDL-3  Dye  Laser)  will  be  used  as  the  uv  source  for  Mie  scattering  imaging  of  the  flame.  The 
Nd:YAG  laser  and  the  pulsed  dye  laser  will  be  explored  as  the  excitation  laser  source  for  LIF 
thermometry.  The  planar  Mie  scattering  of  the  wick  Li-SF6  flame  should  shed  some  light  on  if 
there  exist  condensed-phase  products  in  the  fuel  rich  regime.  The  Nd:YAG  laser  is  also  equipped 
with  a  double  pulse  option  which  can  be  used  to  develop  particle  tracking  velocimetry  to 
characterize  the  hydrodynamic  aspects  of  the  flame.  The  new  diagnostics  will  add  an  additional 
dimension  to  quantify  thermochemical  properties  and  hydrodynamic  characteristics  of  the  wick 
combustion  of  Li  and  SF6,  providing  a  better  understanding  of  the  combustion  process. 
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